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ABSTRACT 

Cosmological simulations of dark matter structures have shown that the equilibrated dark matter 
structures have a fairly small angular momentum. It appears from these N-body simulations that 
the radial profile of the angular momentum has an almost universal behavior, even if the different 
dark matter structures have experienced very different formation and merger histories. We suggest a 
perturbed Jeans equation, which includes a rotational term. This is done under a reasonable assumed 
form of the change in the distribution function. By conjecturing that the (new) subdominant rotation 
term must be proportional to the (old) dominant mass term, we find a clear connection, which is in 
rather good agreement with the results of recent high resolution simulations. We also present a new 
connection between the radial profiles of the angular momentum and the velocity anisotropy, which is 
also in fair agreement with numerical findings. Finally we show how the spin parameter A increases 
as a function of radius. 

Subject headings: galaxies: halos — methods: analytical — dark matter — galaxies: kinematics and 
dynamics — galaxies: general — galaxies: structure 



1. INTRODUCTION 

Our understanding of dark matter structures has in- 
creased significantly over the last years. This progress 
has mainly been driven by pure dark matter numer- 
ical simulations which have suggested or identified a 
range of universalities. One of the first general prop- 
erties to be sugges t ed is the radial d ensity profile 
(iNavarro et alj 119961: IMoore et all 1 19981 iMerritt et all 
120061: IG raham et al.ll2006h . Also more complex connec- 
tions relating integrated quantities have been suggested, 
including the pseudo phase-space dens ity being a power- 
law in radius (jTavlor fc Navarroll2001h . p/a 3 ~ r~ Q , or a 
connection between the velocit y profile and the density 
slope ([Hansen fc Moorel 120061 : iHansen fc Stadell I2006D . 

(3 ~ 7 , where 
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Finally, 



a connection be tween angular mom entum and mass has 
been identified (|Bullock et al.ll200lD . j ~ M s . Only few 
attempts have been made at identifying u niversalities in 
the actual velocity distr ibution function (jHansen et alj 
120061 : IWoitak et al.l l2008h . since non-integrated quanti- 
ties require a very large number of particles in the equi- 
librated structure. 

A wide range of theoretical ideas and models have been 
devised trying to explain these phenomenological pro- 
files and relations. Fi rst of all, the general properties 
of th e density profiles ([Navarro et al.lll996l : IMoore et alJ 
[19981) can be derived analytically under the assump- 
tion that phase-space dens i ty is a power-law in ra- 
dius dTavlor fc Navarr^l2001tiHansenll2004l : lAustin et all 



120051: iDehnen fc McLaughlin! 120051 ). This is done sim- 



ply by inserting the phase-space density into the Jeans 
equation, and then solving it. This procedure reveals 
one unique mathematical solution, which is physically 
plausible in the equilibrated region of the dark matter 
structures. 



A completely different approach is m ad e in 
the Barcelona model (ISalvador-Sole et alj 119981 : 
iManrique et~all 120031 : ISalvador-Sole et alJ 1200*1 . where 
slow accretion in a generalized Press-Schechter model 
allows one to derive density profiles which are in very 
good agreement with the profiles observed in numerical 
simulation. 

It thus appears that there are (at least) two completely 
different possible explanations for the structural prop- 
erty of the density: one being that the density profile is 
slowly grown according to the parameters in the expand- 
ing universe, and the other possibility is that irrespective 
of how the structures are formed, then some unknown 
dynamical process forces the phase-space density to be 
a power-law in radius, which thro ugh the Jeans equa- 
tion gives the density profile itself. lAustin et al.l (2005) 
argue that this process is violent relaxation. A differ- 
ent theoretical approach to get the density profile is a 
series expansion of the coupled collisionless Boltzmann 
and Poisson equ ations which is renormalized in 'time' 
(Hcnrikscn 2007). However, when comparing to numeri- 
cal results, it appears that this last theoretical appr oach 
still loaves room for improvements (Hcnrikscn 2007). 

The relation between angular momentum and 
mass, which was discov ered in numerical simula- 
tions (Bu llock et alj I2001D , has been explored in the 
Barcelona model of accretion-driven formation of cos- 
mological structures (| Gonzalez- Casado et al.l I2007T ) . It 
appears that the results of this theoretical model are in 
rather good agreement with the results of numerical sim- 
ulations. It is therefore tempting to conclude, that the 
kinematic and structural parameters may indeed be gov- 
erned by the physics behind the Barcelona model - at 
least until we find another convincing derivation of the 
j — M relation based on different assumptions. This is 
exactly what we set out to do in this paper. 
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In general when dealing with the angular momentum of 
dark matter in galaxies people often use the so-called spin 

parameter A = > introduced by I Peebles! (|1969t ). as 

a common reference point. Here L, E, M and G are 
the angular momentum, the binding energy and the to- 
tal mass of the system and the Newtonian gravitational 
constant respectively. The spin parameter roughly cor- 
responds to the ratio between the overall angular mo- 
mentum of the object and the angular momentum this 
object needs to sustain rotational support. One of the 
reasons why A is so widely accepted as a good indicator 
of the angular momentum, is that it has the ability of 
being (almost) constant in time, under the assumptions 
that the system is more or less isolated and that there is 
no dissipation present (assuring that E and L are both 
conserved). The us ual size of the spin p arameter is in 
the ballpark of 0.05 (jVitvitska et al.ll2002ft . meaning lit- 
tle systematic rotation and negligible rotational support. 

There is still no agreement on the origin of the angu- 
lar momentum. Two differen t scenarios have bee n tested 
against each other (e.g. by iMaller et all (|2002D ). The 
first scenario states that the angular momentum of struc- 
tures originates from the merging history, and hence is 
dependent on how an d when the given structure formed 
(jVitvitska et al.ll2002D . The second scenario explains the 
appearance of angular momentum as a consequence of 
linear tidal torques between densi ty fluctuations in the 
early stages of galaxy formation (|D'Onghia fc Navarro! 
[2003). 

We will show below that the generalized Jeans equa- 
tion includes a term governing the profile of the angular 
momentum. This allows us to suggest an angular mo- 
mentum profile directly from the Jeans equation, with 
no reference to the specific way the structures were as- 
sembled. This contrasts the claim that merger history is 
crucial when describing angular momentum. We use re- 
sults from recent high resolution N-body simulations to 
show that the angular momentum profile of the numer- 
ically equilibrated dark matter (DM) structures shows 
good agreement with our suggested relation. This indi- 
cates that the angular momentum profile is fixed through 
the Jeans equation, irrespective of the detailed struc- 
ture formation proce ss. This result is suppo r ted b y the 
numerical results of lAscasibar fc Gottldberl (|2008l ) and 
iD'Onghia fc Navarro! (|2007ft . who conclude that "equilib- 
rium dark matter haloes show no significant correlation 
between spin and merging history" . 

We present a connection between the radial profiles 
of the angular momentum and the velocity anisotropy. 
This correlation is also shown to be in fair agreement 
with numerical simulations. 

Finally we combine our suggested relation between 
the angular momentum and the mass with the spin pa- 
rameter, and find that the radial dependence of this 
spin parameter is in fair agr eement with the results of 
lAscasibar fc Gottldberl (|2008l ). 

2. JEANS EQUATION INCLUDING ROTATION 

In general nothing ensures that the particle ensemble 
in a DM halo is spherical. Nevertheless it makes the 
approach relatively simple and analytically manageable, 
without being far from the 'real' triaxial picture. The 
system is then governed by the spherical(ly symmetric) 



collisionless Boltzmann equation (CBE) (with the veloc- 
ity distribution function being a function of the radial 
dependence alone). One might a lso use other (more com- 
plicated) forms of the DF (e.g. iTonini et al.l (|2006l )) but 
this would make the approach analytically cumbersome. 
The CBE describes the relation between the spatial dis- 
tribution function, the gravitational potential and the 
velocity distribution of the DM particles. 

Combining the first moment of the CBE with the 
spherical Poisson equation, one can write the general 
spherical Jeans equation for the DM particle ensemble, 
under the assumption that the thermal velocities are in- 
dependent of eac h other, and that there is no bulk (ro- 
tational) motion (|Binnev fc Tremaindll987l ) 
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Here Og, of and are the velocity dispersions, p is the 
density, M is the total mass within radius r, and 67 is 
the gravitational constant. 

We wish to consider the angular momentum, so we 
will now include a small radial dependent bulk rotation, 
v \:ot( r )- W e explicitly write the velocity in the following 
way 

v 4>( v mt = 0) -> ^K-ot) . ( 2 ) 

which means that the azimuthal velocity goes from be- 
ing purely thermal to containing both thermal and bulk 
motion. 

In general v TO i is not a function of the spherical radius 
r, but more likely a function of the cylindrical radius 
R. So when we in the following refer to v r0 ^(r), what we 
really mean is, that the radial velocity has been averaged 
over a spherical shell with radius r. This is also the 
way v ro t is defined in the simulations we will compare 
with in section [3] Furthermore since the effect of bulk 
rotation cannot depend on direction, we only consider 
the absolute value of v ro ^. 

The simplest way to add rotation to our system, is to 
give every particle an initial kick. This would result in 
a shift of the distribution function (DF) /, towards a 
higher mean velocity by the amount added, call it w ro t- 
Since we need the structure to be a relaxed system in 
equilibrium we can't allow this. A pure shift of the DF 
would enable the most energetic particles to escape the 
system, and thereby bring it out of equilibrium. Instead 
of just shifting the DF we force it to be distorted into 
a new DF /, making sure that the system is always in 
equilibrium. 

To be able to implement such a distortion into Eq. |T]) 
we will now parametrize the distortion of the DF. As 
a first step, we assume that the difference in the two 
DFs is just a small perturbation to the overall system. 
Furthermore the distortion must depend on the added 
bulk rotation so that we have 



/ - / cx P(v rot )f 



(3) 



Here P is some unknown function depending on v r0 ^ 
which we will evaluate later. 

As mentioned above, simply shifting the DF by v T0 ^ 
while leaving the shape of the DF virtually unchanged 
would allow energetic particles to escape. The pertur- 
bation of the DF must therefore also depend on the az- 
imuthal velocity itself, such that the amount of distortion 
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is not the same at all velocities. Introducing another un- 
known function Q we then have 

/-/ K P( %t )QW/ (4) 

Furthermore, we will demand that the overall density 
of the system is not affected by the added bulk rotation, 
i.e, p = p. The density is statistically defined as 

P = J M 3 v = J /-(/- f)d 3 v = p- J(f- f)d 3 v (5) 

where we have used that / = / — (/ — /)• Ensuring 
that the density is unaffected by the added rotation is 
easily done by restricting (/ — /) to be an odd function 
when integrated over the velocities, such that the last 
integral vanishes. Assuming that Q is a simple power 
law in the azimuthal velocity with a positive odd integer 
exponent (to be motivated later in this section) this is 
accomplished and we have 



/ - / = £P(v. 



rot J 



K) 1/2 



(6) 



with 7 = 2n+l for n = 0, 1, 2, ... and £ being an unknown 
constant. 

The Q(v<p) part of the parameterization on this power 
law form will in principle give problems for v^, — ► — oo, 
causing the distortion of the DF to become infinite and 
hence / to become negative, but since f(v^,) is nearly zero 

for W0 ~ ±4^/cr^, this is not a practical problem. How- 
ever, for a more formal derivation one naturally cannot 
allow / to become negative. 

Thus, by making 7 an odd integer, what we have done 
(using Eq. ((1])) is to create a new DF, /, that generates 
the density p by adding an odd-powered DF concerning 
to our origin al DF, /. This is somewh at similar to 
the discussion in lBinnev fe; Trem ainc (1983) section 4.5. 

In principle, the assumption in Eq. © can be tested 
with high resolution N-body simulations, by consider- 
ing the two tangential velocity distribution functions 
(VDFs), and simply looking for the difference be- 
tween the shapes of the VDFs in the directions par- 
allel and perpendicular to the angular momentum vec- 
tor. This we can do using the N-bo dy/gas dynami- 
cal s i mulated large disc galaxy 'K15 ' (|Sommer-Larsenl 
[2fMlHansen fc Sommer-Larsenl[2006l) . The 'K15' simu- 
lation is a significantly improved version of the TreeSPH 
code used previously for ga laxy formation simulations 
(jSommer-Larsen et al.l [2003^ . This simulation is based 
on a flat ACDM model with (Om.Aa) = (0.3,0.7). The 
simulation takes many of the important factors of galaxy 
formation into account such as SN feedback, gas recycling 
(tracing 10 elements), atomic radiative cooling, etc. It 
consist of both cold and warm gas, DM, disk and bulge 
stars and stellar satellites. The galaxy 'KI5' contains 
about 3 x I0 5 gas and DM particles and it has mgas = 
m stars = 7 - 3 x 10 5 M Q //i and m DM = 4.2 x 10 6 M Q /h 
where h = 0.65. Furthermore the gravitational (spline) 
softening lengths adopted are egas = e stars = ^80 and 

e DM = 680 P c / h - 

In Fig. Q] we have plotted the VDFs of the DM parti- 
cles of one of the bins outside the stellar disk region in 
'K15'. We clearly see that the DF including rotation (red 




Fig. 1. — The velocity DFs for the DM particles in one of the 
potential bins in the galaxy K15 in the radial direction (black solid 
line) and the tangential direction perpendicular to (red dashed line) 
and parallel to (blue dot-dashed line) the angular momentum vec- 
tor for the considered bin. We clearly see that the DF containing 
rotation (perpendicular to the angular momentum vector) is dis- 
torted as suggested. All the trustable bins in 'K15' show a similar 
picture. This supports our assumption of the added rotation caus- 
ing a distortion of the DF of the DM structure. 

dashed line) is indeed distorted compared to the radial 
VDF (black solid line) and the tangential VDF parallel 
to the angular momentum vector (blue dot-dashed line) . 
This supports our assumption of a distortion of the DF 
caused by the added bulk rotation. Furthermore we real- 
ize that the distortion seems to grow (i.e., the difference 
between the two red dashed lines enhances) as the veloc- 
ity grows out to about 2-3 a, beyond which the number 
of particles in each bin of the simulation is too small to 
conclude anything. This further motivates a form of Q 
similar to the one suggested in Eq.©. 

To be able to see the effects the added rotation has on 
the system in general, the goal is to get an expression for 
the Jeans equation concerning a system to which a small 
bulk rotation has been added. This can be accomplished 
by making an expression for the new perturbed velocity 
dispersion in the azimuthal direction, 5^. The velocity 
dispersion of a system without rotation is given by 



fv, 



2 d 3 v 



(7) 



By definition the velocity dispersion is the integral over 
the DF multiplied with the difference between the indi- 
vidual particle velocities and the mean velocity of the 
system. Since the mean velocity after adding rotation 
becomes equal to the added rotation itself, the new per- 
turbed azimuthal velocity dispersion must take the form 



- / fir,, ■ <y,,j T ,t'c 



(8) 



where it is easily shown that v T0 ^ = 1/p J fv^v. 

Using that J fd 3 v — p, i.e., / is an even function in 
velocity space, that (/ — /) is an odd function to conserve 
the density and that /=/—(/ — /), combining Eqs. (0) 
and (H|) gives 



Tot 



fd 3 v- 
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Since it is known from numerical cosmological simula- 
tions that the rotational ener gy is less than a few per- 
cent of the thermal energy ( Bullock et al.l [20011, i.e., 
v rot ^ a % ^ * s justified to ignore the higher order term 
in i> ro t- This combined with Eq. j7]) implies that 



2w rot 



(/ - f)v 4 d 3 v 



(10) 



Thus the new velocity dispersion can be written as the 
old one plus a term concerning the distortion of the DF 
as well as the rotation. Combining this with the DF 
distortion in Eq. ([6]) implies 

al = al + 2Zv I0t P(v T0t )(al)-l~ J vfa f d 3 v . (11) 

We are then left with evaluating an integral on the form 
1 



P 



(12) 



where k — 7 + 1 is an even number (gamma is odd). 

Recognizing that the integral in Eq. (fT2"|) is just the ex- 
pression for the fc'th moment, a^, for v^, centered around 
the mean (v^) = gives 

- J v$f d 3 v = a k {al) k/2 with k = 2n+2 for n = 0, 1, 2, .. 

Combining this with the expression for the perturbed 
velocity dispersion in Eq. (|llj) gives 



^•2 _2 

aj, = a j, 



2£v Iot a 1+ i?(» rot ) 



(14) 



We are now able to quantify the v IQ ^ dependency of 
the distortion of the DF, i.e., the function P. An easy 
way to evaluate the function P is by rearranging Eq. JTj 
so that 



p (^rot) 



(15) 
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Using the 'K15' simulation again we are then able to es- 
timate the actual size of P as a function of the rotation. 
Calculating the relevant quantities from the simulation 
and plotting P as a function of v T0 ^ gives Fig. [2] From 
this figure, which is showing the numerically resolved re- 
gion of the structure, we see that P seems to depend 
liniarly on the rotation, meaning that in this region of 
the structure P — C% + C*2V ro t where the Cs are con- 
stants. However one must keep in mind that (outside 
the resolved region) P must go to for vanishing v I0 ^. 
Combining this with expression (| 14[) and again ignoring 
higher order terms in v TO i leaves us with 



°& = °l + V 12 a 7 +i «rot 



(16) 



Here we have introduced the constant r\ = (for sim- 
plicity). Thus, what we have done here is simply insert- 
ing the linear dependence of P on v ro i into Eq. (fT4|) and 
using that ^ <C ct| as done in Eq. (fT0|) . 

As mentioned a is just the moment corresponding to 
the chosen value of 7 (where Q ~ v^), i.e., a constant. 
And since the above result only relies on the restrictions 
on 7 it implies that the choice of 7 doesn't result in loss 



Fig. 2. — The function P(v r nt) from the simul ation 'K15' 
(Sommcr-Larscn 2006; Hansen & Sommer-Larsen 2006) plotted as 
a function of v ro ^ in the trustable region, i.e., from around 10 kpc 
(where most of the gas disk vanishes) and out to the virial radius. 
For f ro t — > one must have P — > which however is outside the 
numerically resolved region in this specific case. 

of generality. Thus we are free to choose any value of 
7 when investigating the above equation. We will chose 
7 = 3 since we are then able to estimate the size of 
the moment. The fourth moment, i.e., the kurtosis of a 
Gaussian DF is 3, and since we expect the DFs of DM 
structures to be Gaussian like, using 7 = 3 will make it 
easier to compare with simulations. This implies that 



al +?7l2 a 4 w ro t \ 



(17) 



This expression is of course a consequence of the assumed 
form of Q. One could definitely argue for other forms of 
Q fulfilling the request p = p, e.g. an exponential form 
or a combination of both exponential and power law. 
However in order to simplify the analytical calculation of 
the integral in Eq. (|10[) we have chosen the simple power 
law form. In the future it would definitely be interesting 
to test other forms of Q (and P) to see if this effects 
the final results and conclusions significantly. We intend 
to make a qualitative estimation of both Q and P in a 
following paper. This will, among other things, require 
a larger sample of high resolution equilibrated structures 
(both pure DM as well as DM+gas simulations) and that 
we systematically test the effects of using either spherical 
or potential bins, of the structures non-sphericity etc. 

Combining Eq. (| 1 T|) with Eq. ([1]) (for a perturbed sys- 
tem, i.e., — ► al) leaves us with a Jeans equation con- 
taining four terms. The 3 normal ones, dealing with the 
density, mass and velocity dispersion profile, and one new 
term describing the effect that rotation has on the system 
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where the anisotropy is given by = 1 — ag/a 2 . Here 
we assume that a# = cr?, which states that the thermal 
velocity moments in the tangential directions are equal, 
irrespective of the magnitude of the (small) bulk rota- 
tion. Note that simply substituting cr^ into the (per- 
turbed) Jeans equation leaves us with a Jeans equation, 
only involving the unperturbed velocity dispersions. 

The dominating terms in Eq. fTS]) are the derivative 
and the mass terms. Making the conjecture that the 
rotation term, which is just a minor perturbation of the 
Jeans equation, must follow the profile of the dominating 
mass term, and assuming (for now) that /? = 0, we get 
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directly from Eq. (| 18[) a relation between the rotational 
perturbation and the dominating mass, which reads 



GM 



v ToV 



(19) 



In principle many other solutions, than the conjecture of 
the small term following the dominant one used above, 
are allowed to exist, but these would all imply some de- 
gree of compensation or fine-tuning between the various 
terms. We therefore suspect that there is a more physical 

explanation for why the v T0 ^ \Jo^, term is proportional to 

than our conjecture, but none has been found so far. 

The different structures may have fairly different mag- 
nitudes of the angular momentum, and Eq. (|19| expresses 
only that the radial profile of the angular momentum is 
always the same, however, the absolute magnitude is un- 
known, and may vary from structure to structure. 

In a similar way we can look at the relationship be- 
tween the (3 ^ and the rotational term. We find from 
Eq. (TT8| that this connection is 



(20) 



This relation implies that if (3 goes to 0, the rotation 
term should go to as well. Since we are here suggesting 
a relation between the two minor terms in Eq. (|18|) the 
relation (|2H)l might not be as strong as relation (TT§|) . The 
case /3 < does not occur in the equilibrated part of 
the simulated DM halo structure and has therefore no 
relevance to the problem at hand. 

We are aware that p is marginally smaller than in 
DM05. In fact some of our structures also have f3 < 
in some of the inner most bins. But since we, as well as 
DM05, are working with simulations which are known to 
have difficulties simulating structures at the innermost 
parts, such values (which are not much below 0) must be 
considered in agreement with within errors and does 
therefore not conflict our suggested relation between beta 
and u ro t- If on the other hand cosmological simulations 
were to produce an equilibrated structure with a clear 
trend that a fully resolved smooth region of the structure 
had (3 < this would definitely question the validity of 
our work. 

Note that including a centrifugal term into the equa- 
tions will basically give a small energy conserving per- 
turbation, which goes as "^of ^° ^ e new azinlu th a l vc ~ 
locity dispersion. However this perturbation is so small 
that it is not visible in the numerical simulations, and it 
is therefore ignored. 

It is now straight forward to test these suggested rela- 
tions with the results from numerical simulations. 

3. COMPARING WITH NUMERICAL SIMULATIONS 

We have argued that there may be clear relations be- 
tween the new rotational supplement to the Jeans equa- 
tion and the mass- and anisotropy-terms. To test this 
we used 10 intermediately resolved galax y and cluster 
sized numerical simulations of DM halos (|Macci6 et al.l 
12001 . one high resolution cluster, C#a..W3, and one 
high resolution galax y, the 'Via Lactea' simulation 
(jDiemand et al.ll2007al fbl). 

The 10 intermediately resolved simulations have been 
performed using PKDGRAV, a treecode written by 




4 6 
log( 4.30*1CT 5 *M(<r) ) 

Fig. 3. — The suggested connection between the mass (r.h.s. 
of the Jeans equation) on the x-axis, and the new rotational term 
(solid line) on the y-axis. The diamonds, triangles, crosses and the 
squares represent the ga laxy sized halos and the cluster sized ha- 
los from the simulations Maccio et al. (2007), the high resolution 
C£ffl.W3 cluster and the Dicmand ct al. ( 2007a) "Via Lactea" high 
resolution simulation respectively. This illustrates that the Jeans 
equation determines the radial behavior of the rotation, i.e. the 
angular momentum j(r). We use SI units and 04 = 3. The factor 
4.30 ■ 10 -6 includes the gravitational constant G, and is the cor- 
rection needed to have the quantities in SI units. The rj for each 
structure corresponds to the ones listed in Table [T] 



Joachim Stadel and Thomas Quinn (|Stadell l2001). The 
initia l conditions are ge nerated with the GRAFIC2 pack- 
age ([Bertsc hingcr 2001). The starting redshifts Zi are set 
to the time when the standard deviation of the smallest 
density fluctuations resolved within the simulation box 
reaches 0.2 (the smallest scale resolved within the initial 
conditions is defined as twice the intra-particlc distance) . 
All the halos were i dentified using a SO (Spherical Over- 
density) algorithm (iMaccio et al.ll2007h . The cluster- like 
haloes have been extracted from a 63.9 Mpc/h simula- 
tion containing 600 3 particles, with a mass resolution of 
m p — 8.98 x 10 7 Mq/Ii. The masses of the clusters used 
for this study are 2.1, 1.8, and 1.6 xlO 14 M //i. The 
galaxy sized haloes have been obtained by re-simulating 
at high resolution haloes found in the previous simu- 
lation. The simulated haloes are in the mass range 
0.9 — 2.5 x 10 12 Mq//i and have a mass resolution of 
m p = 4.16 x 10 5 Af Q //i that gives a minimum number of 
particles per halo of about 2.5 x 10 6 particles. The High 
resolution cluster Ch_rW3, based on the PKDGRAV as 
well, has 11 millions particles within its virial radius and 
a mass of M = 1.81 x 10 14 M Q /h. The 'Via Lactea' 
(which is also based on the PKDGRAV code) simulation 
includes 234 million particles with force resolution of 90 
pc, and it includes one highly equilibrated structure of 
mass M2o n = 1.77 x 10 12 MfT)i con taining about 84 million 
particles (|Diemand et al.ll2007a| ). 

Plotting the rotation-term against the mass for the dif- 
ferent simulations gives Fig. [3) Here the diamonds, tri- 
angles, crosses and squares represent the galaxy sized 
halos, t he cluster sized halos, the C#_r.W3 simulation 
and the iDiemand et all (|2007a[ ) "Via Lactea" high res- 
olution simulation respectively. In Fig. [3] we see a clear 
linear relation between the two terms. This means that 
the generalized Jeans equation (Eq. (fT8|) ) determines the 
radial behavior of the rotation, i.e. the angular momen - 
tum j(r). This also explains why iBullock et all (|2001f ) 
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Fig. 4. — The relation between the thermal velocity anisotropy /3 
and the rotation-term (straight lin e). The squares are from the high 
resolution numerical simulation of Dicmand et al. (2007a) and the 
diamonds are the galaxy sized halos from the Maccio ct al. (2007) 
simulations. The crosses is the Cjj\r.W3 high reso lution cluster 
and th e triangles are the cluster sized structures from Maccio ct al. 
(2007). We see an indication of a deviation from the relation for 
small j3. This is probably a consequence of comparing the two 
subdominant terms in the new Jeans equation. When plotting 
the data we have cut off the structures so that all of the points 
can be considered to be equilibrated. We have determined the 
cutoff-points by excluding the (outer) part of the structures, where 
the velocity anisotropy is not a (roughly) monotonically increasing 
function of radius. Furthermore we have re-binned the data to 
reduce scatter, so that each point now contains roughly 1/7 of the 
structure. We again use SI units, 0:4 = 3 and different r\ for each 
structure as listed in Table [T] 

find a strong relation between the angular momentum 
and the mass in their simulations, since ou r conjecture 
resembles the results of iBullock et al.l (|2001h when er^ is 
constant. We have tested that this relation is not just 
an effect of choosing (actually deriving) a term in the 
Jeans equation with the right units. For instance the 
term wit h t^nt 7 " does not have a correct relation to the 
mass (aslHoisgaar d et all (|2007l ) also conclude). 

In a similar way we can test our suggested linear re- 
lation between of/3 and f r oty^|- Plotting these quan- 
tities for the intermediate resolution halos together with 
the C/ffl.W3 and 'Via Lactea' high resolution simula- 
tions gives Fig. |U Here we see a clear correlation for 
the majority of (3 values. However, there is some indi- 
cation that small (3 values doesn't follow our relation as 
strictly as for large (3 values. This is probably a conse- 
quence of comparing the two subdominant terms in the 
new Jeans equation with one another, which as men- 
tioned doesn't make relation (|20|) as strong as the re- 
lation between the mass and the bulk rotation. In the 
figure we have re-binned the data to reduce scatter and 
cut off the structures where (3 was no longer a (roughly) 
monotonically increasing function of radius. Plotting the 
structures without making an outer cut doesn't change 
the picture but only enhances the overall scatter. 

We see that the rotation term goes to as (3 goes to 
0, exactly as suggested. In fact, if we plot the fraction 
of the kinetic energy in rotation, i.e., wj^/cr^, we see 
that it drops from 10~ 2 in the outskirts of the structure, 
down below 1CP 4 for the inner-most bins. Since (3 is 
monotonically increasing as a function of radius, the fact 
that the rotation becomes so small in the inner parts 



Fig. 5. — The free parameter »7mass obtained from the relation 
plotted in Fig.[3]and rip obtained from the relation in Fig. [4] plotted 
against each other. The solid guide-the-eye line represents the 
relation »7mass = 2.577,3. The error bars represent the interval in 
which the r\ values give the best fit to the relation. For instance 
the 'Via Lactea' (square) fits the mass relation in Fig. [^reasonably 
well for 0.17 < rj mass < 0.47, i.e., 0.32±0.15 as written in Table.[T| 
The symbols are the same as in Fig. [3] 

of the structure supports our suggested relation of the 
rotation term going towards for small (3. 

In Figs. [3] and 0] the only free parameter in our rela- 
tions, 77, have been fitted for each structure. These val- 
ues of 77 corresponding to the relations in Eqs. (fl9|) and 
(|20p represents the unknown magnitude of the angular 
momentum and are shown together with the estimated 
errors in Table [TJ 

Plotting the 77 values and their errors gives Fig. [5] Here 
we see a tendency of ?7 m ass being larger than 77^3 . Since 
the triangles and the cross are cluster like structures and 
the rest are galaxy like structures we notice that there 
might be a connection between 77 and the mass of the 
structures. In Fig. [6] we plot rjp (since it has the smallest 
error bars, percentage-wise) against the estimated virial 
mass of each structure (see Table [TJ , and see that our 
free parameter anti-correlates slightly with the mass of 
the structure. So according to our suggested relations 
the effect an added bulk rotation has on a system is anti- 
correlated with the virial mass of that system. 

After having tested our suggested relat ions from the 
previo us section with the simulations from iMaccio et ail 
(|2007f ). we also held them up against the recent 
high resolution numer ical simulation 'Via Lactea' by 
iDiemand et al.l (|2007af ). to make sure that the results 
is not just a coincidence because of lack of numerical 
resolution. We have plotted the high resolution data as 
squares in Figs. [3] to [6] for comparison. In Figs. [7| and 
[5] we have plotted the 'Via Lactea' alone, without any 
cutoffs or re-binning. In these figures we are using the 
values 77 from Table [TJ On both figures we see that the 
suggested relations are confirmed when comparing with 
highly resolved data. The 'Via Lactea' structure did not 
experience any major mergers since z = 1. All quanti- 
ties are extracted in spherical bins. Due to numerical 
softening one can safely trust the radius outside 1 kpc of 
this galaxy. In this simulation the outermost 6-10 points 
should be considered with care since they are potentially 
not fully equilibrated yet, as is easily seen when consid- 
ering the radial derivative of the density profile. 
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TABLE 1 Our free parameter 17 from Figs. [3| and [4) its errors and the mass 

AND AZIMUTHAL VELOCITY DISPERSION TAKEN AT THE VIRIAL RADIUS r vir , FOR THE 



STRUCTURES USED TO TEST OUR RELATIONS. THE ERRORS ON T\p AND r) maS s 
REPRESENT THE INTERVAL IN WHICH THE VALUES GIVE THE BEST FIT TO THE 
RELATIONS. FOR INSTANCE THE 'VlA LACTEA' FITS THE MASS RELATION IN FlG. QJ] 
REASONABLY WELL FOR 0.17 < ?Jmass < 0.47, I.E., 0.32 ± 0.15. ALL THE STRUCTURES 

ARE BASED ON THE PKDGRAV TREE CODE. THE 'VlA LACTEA' AND THE GX.XX 
ARE GALAXY SIZED STRUCTURES AND THE C_ff_R.W3 AND CX.XX ARE CLUSTER SIZED 
STRUCTURES. THE XX. Wl IS BASED ON A WMAP 1-YEAR DATA COSMOLOGY 

(h = 0.71 and Q.m = 0.268), whereas XX. W3 is based on a WMAP 3-year data 
COSMOLOGY (h = 0.73 AND Q m = 0.238). For more information o n the 
structures, refer to the text or se e Diemand et al.I (|2007aD and 
IMaccio et al.I (120071) . 
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Fig. 6. — The free parameter r\p plotted against the virial mass 
of each structure (listed in Table [T] The symbols are the same as 
in Fig. [3] 

We have thus compared the angular momentum of 
the intermediate resolution structures, the highly re- 
solved C.ff_R.W3 cluster, and the 'Via Lactea' simulation, 
whic h is one of the best reso lved structures published to- 
day (|Diemand et al.ll2007a[ ). to our suggested relations, 
and see strong correlations between the rotation, mass 
and velocity anisotropy of the system. 

As mentioned in the introduction, people often use the 
spin parameter A when describing the angular momen- 
tum of DM halos. Combining our relation between the 
mass and the ang ular momentum with the spin parame- 
ter, as defined bv lBullock et alj (|2001[ ) 



A' = 



J 



(21) 



V2GM 3 r 

we end up with a new expression for the spin parameter 

1 Vr 



A' = 



V2 j] 12 a 4 



(22) 



Fi g. 7. — The same plot as Fig. \3\ containing only the data from 
the Diemand et al. (2007a) "Via Lactea" high resolution simula- 
tion . We only plot numerically simulated points (squares) in the 
resolved region, i.e. outside r = 1 kpc, and out to r = T200- The 
outermost 6 to 10 points (top right corner) are possibly not yet 
fully equilibrated (as can be seen when looking at the profile of 
an d might therefore be ignored. We use SI units, r\ = 0.32 
and 0:4 = 3. The factor 4.30 ■ 10 — 6 includes the gravitational con- 
stant G, and is the correction needed to have the quantities in SI 
units. 



the spin parameter only as a function of mass and 

without any dependence on the bulk rotation i> ro t- In 
Fig. [5]we have plotted Eq. (f2"2"| for the simulated struc- 
tures. This figure agrees f airly well with Fig. 5 of 
lAscasibar fc Gottloberl (|2008f ). If we estimate a gradient 
of our plot we get approximately 1/4 to 1/5 (depending 
on the chosen s tructure), whereas an estimat ed gradi- 
ent on Fig. 5 in Ascas ibar fc Gottloberl (|2008l ) is closer 
to 1/6. Thus Eq. (fi"9]) appears to roughly explain the 
observed tendency of an increase in the spin parameter 
as a function of radius. As mentioned earlier a relation 
with a rv^ ^ term instead of the rv r0 ^ J a\ term we sug- 



where v 2 c = GM/r. We are therefore able to describe gest will according to simulations not give the correct 
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Fig. 8. — The relation between the thermal velocity anisotropy 
(i and the rotation-te rm (straight line), and t he high resolution nu- 
merical simulation of Dicmand et al. (2007a) (squares). The trend 
that lower j3 doesn't obey our relation as well as large j3 is also 
obvious here. We again use SI units, r] = 0.26 and 04 = 3. The 
outermost 6 to 10 points are possibly not yet fully equilibrated 
(clearly visible when considering as a function of radius) and 

might therefore be ignored. 
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Fig. 9. — The spin parameter (Eq. i22l ) plotted as a function of 
radius. The symbols are the same as in Fig. [3] 

relation to the mass. Furthermore a result on the fv^ Q ^ 
form implies a constant spin parameter as a function of 
radius, and since this is clearly no t in ag reement with 
the work by lAscasibar fc Gottfoberl (|2008l ). we take this 
as yet another indication of the success of the relations 
in Eqs. (fl9 |) and ([20]) . 

Finally it have been suggested that the spin parame- 
ter doesn't depend on the virial mass of the structures 
(|Macci6 et al.ll2007l ). To test this we plot in Fig. [TUl 
the values of A' from Eq. (f2"2"|) taken at the virial radius 
of the structures. To calculate A^- r we used the values 
given in Table [TJ Here we see the indication of a slight 
increase in spin as the virial mass of the structures grow. 
The linear fit in Fig. \W\ (full line) has an inclination of 
0.09 ± 0.04. One possible explanation for the indication 
of a mass dependence, might be the fact that the spin is 
taken at r v ; r . As pointed out by lAscasibar fc Gottloberl 
(2008) the use of r vir (compared to their i? max ) might be 
too 'non-conservative' when estimating the various prop- 
erties of equilibrated structures. On the other hand it is 
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Fig. 10. — The spin parameter (Eq. i22l ) at the virial radius. 
The linear fit to the points (straight line) is given by (0.09 ±0.04) X 
log(Af v ; r /M0) — (2.31 ± 0.51). The errors correspond to the error 
in r;mass (see table [TJ used when calculating X' vir - The dotted lines 
repr esents the mean (01 )34 ± 0.001) and 2a scatter (0.55 ± 0.01) 
from Maccio et al. (2007), which are in agreement with the results 
of Bullock et al. (2001). The symbols are the same as in Fig. [3] 

not surprising with a slight increase in spin, since the 
spin parameter as defined in Eq. I122|) basically resem- 
bles a relation between ?7mass and -Myir Sirn il ar to the 
one shown for 77^ in Fig. [5J Nevertheless, because of the 
large scatter and error-bars in our points we must con- 
clude that there is no (significant) dependence between 
our spin parameter and the virial mass of the structures. 
This is in agreeme nt with the lower part of Fig. 3 in 
iMaccio et al.l (|2007|) , which also shows that a relatively 
large scatter in the points is usual. 

4. CONCLUSIONS 

We have studied the form of the spherical Jeans equa- 
tion when one includes angular momentum, and we find 
that a bulk motion leads to the introduction of an ex- 
tra term, which includes the average rotational velocity. 
This is done under the assumption that the distortion 
of the distribution function takes the form argued for in 
Eq. @. This assumption is supported by numerically 
simulated structures. The distortion enables us to sug- 
gest a new correlation between the angular momentum 
and the mass. This relation is in good agreement with 
the findings of recent high resolution numerical simula- 
tions of cosmological structures. We also suggest a corre- 
lation between the angular momentum and the velocity 
anisotropy profiles, which is also in fair agreement with 
numerical findings. These suggested relations imply that 
cosmological dark matter structures have angular mo- 
mentum profiles which have the same universal proper- 
ties, irrespective of how or when they were formed. 

Finally we derive a new form of the spin parameter, 
A', which is shown to increase slowly as a function of 
radius, in agreement with recent simulations. Further- 
more our relations indicate that there is no (significant) 
dependence between X' vh . and M v : r . 
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